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Abstract. 

We study the two-spin entanglement distribution along the infinite S — 1/2 chain 
described by the XY model in a transverse field; closed analytical expressions are 
derived for the one-tangle and the concurrences C r , r being the distance between 
the two possibly entangled spins, for values of the Hamiltonian parameters close to 
those corresponding to factorized ground states. The total amount of entanglement, 
the fraction of such entanglement which is stored in pairwise entanglement, and the 
way such fraction distributes along the chain is discussed, with attention focused on 
the dependence on the anisotropy of the exchange interaction. Near factorization 
a characteristic length-scale naturally emerges in the system, which is specifically 
related with entanglement properties and diverges at the critical point of the fully 
isotropic model. In general, we find that anisotropy rule a complex behavior of the 
entanglement properties, which results in the fact that more isotropic models, despite 
being characterized by a larger amount of total entanglement, present a smaller fraction 
of pairwise entanglement: the latter, in turn, is more evenly distributed along the chain, 
to the extent that, in the fully isotropic model at the critical field, the concurrences 
do not depend on r. 
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1. Introduction 

The analysis of entanglement properties has recently furnished new insights into several 
peculiar features of many-body systems, such as the occurrence of quantum phase 
transitions, or that of non trivial factorized ground states [1-11]. Different types of 
entanglement can be defined in many-body systems, but computable measures are 
available just for a few of them. In this sense, a privileged role is played by the 
bipartite entanglement of formation [12], and by the related quantities, one-tangle and 
concurrence, which represent the entanglement of formation between one qubit and the 
rest of the system, and that between two selected qubits of the system, respectively. In 
particular the definition of the concurrence holds not only for pure states, as in the case 
of the one-tangle, but also for mixed ones [13, 14]. 

When magnetic systems are considered, the qubit is naturally represented by a 
spin with S = 1/2: for interacting magnetic models described by Hamiltonians with 
certain symmetry properties, both the one-tangle and the concurrence are expressed in 
terms of standard magnetic observables, such as the magnetizations and the correlation 
functions, making it feasible a quantitative analysis of the entanglement dependence on 
the Hamiltonian parameters. If analytical expressions are available, a general discussion 
of such dependence is at hand, which is the reason why a renewed interest is being 
devoted to exactly solvable models. 

In this paper we focus on the T = behavior of the best known one- dimensional 
XY model in a transverse field, in the vicinity of factorized ground states [10,15]. We 
derive closed analytical formulas for the magnetization and the correlators, as functions 
of the anisotropy, the field, and the distance between the two selected spins, which 
allow us to study the long-distance behavior of the concurrence as the Hamiltonian 
parameters are varied. By using these expressions for the correlation functions, we 
explicitly proof the divergence of the range of the concurrence in the anisotropic model 
(whose preliminary result was presented in Ref. [16]) and we extend it in the whole 
parameter space, studying the fully isotropic case as well as the slightly anisotropic 
region close to the spin saturation. Our analysis, besides analytically confirming the 
divergence of the range of the concurrence, shows that to such divergence corresponds 
the appearance of a characteristic length-scale in the system, that we have named two- 
spin entanglement length. This length-scale depends on the value of the anisotropy and 
keeps finite as far as the model belongs to the Ising universality class, while diverging 
when factorization gets to coincide with saturation, i.e. for the isotropic XX model, 
which belongs to the Kosterlitz-Thouless universality class. Correspondingly, the way 
the two-spin entanglement distributes along the chain is found to strongly depend on 
the symmetry of the model: in the XX model a good amount of entanglement can be 
stored even between two spins which are far apart from each other, while in the Ising 
model the pairwise entanglement of the ground state is shared only between nearest and 
next-nearest neighboring spins. 

The structure of the paper is as follows: In Sec. [2] we introduce the model and 



Two-spin entanglement distribution near factorized states 



3 



the entanglement properties we aim at studying; In Sec. [3] we study the long-distance 
pairwise entanglement both in the anisotropic and in the isotropic case; In Sec. H] we 
define and analyze the two-spin entanglement length, while in Sec. [5]we use our results to 
understand the interplay between pairwise entanglement and multipartite entanglement. 
Finally, in Sec. El we draw the conclusions. 

2. Model 

The XY model in a transverse field is described by the Hamiltonian 

H = J H K 1 + 7)5f5*.i + (1 - l)SfS! +1 - hSf] , (1) 

i 

where i runs over the sites of an infinite chain, Sf (rj—x, y, z) are the 5 = 1/2 quantum 
spin operators, 76 [0,1] is the anisotropy, and h = g^sH/J is the reduced magnetic 
field; J > is the strength of the exchange interaction. 

For < 7 < 1 the model belongs to the Ising universality class and at T = 
the critical field h c — 1 separates a disordered phase (h > h c ), from a spontaneously 
broken-symmetry phase, where the staggered order parameter is finite ((Sf) 7^ 0). In 
the isotropic case, 7 = 0, the model has an addictional rotational symmetry on the 
xy plane, and the critical field coincides with the saturation field, above which all the 
spins incoherently align parallel to the field. For 7 = and h > 1 the system is in a 
fully-polarized phase ((Sf) = |), while for h < h c , the systems is in a gapless phase 
with (Sf) < |, (Sf) = and power-law decaying correlation functions in the xy plane. 
No spontaneous symmetry breaking is present in the isotropic testified by (Sf) 

being null for whatever value of the applied field. 

Let us now consider the h — 7 parameter space of the Hamiltonian ([1]) : The ground 
state of the model is exactly factorized [15] 

\GS) = l[\<p i ) , (2) 

i 

on the circle h? + 7 2 = 1, as well as along the line {h > 1,7 = 0}: Such ground 
state has a Neel structure given by = (— l^cos&yl fi) + sin# 7 | If) , with 

cos# 7 = a/(1 — 7)/(l + 7) = a, which reduces to the trivial ferromagnetic ground 
state for 7 = and h > 1. In what follows, we will refer to the circle h 2 + 7 2 = 1 as the 
factorized circle, and to the line {7 = 0, h > 1} as the factorized line. 

For the model Eq. (jTJ) the concurrence C r between two spins sitting on sites i and 



j, with |i — j\ = r, reads [17] 

C r = 2max{0,C;,C;}, (3) 

C' r = \gf x + gf\ - \j(\ + 9r?-Ml (4) 

C':=\gr-9 y r y \+9 Z r Z -\, (5) 

while the one-tangle is T\ = 1 — 4(M 2 + Mf) , with the correlators = (S?Sj) and 



the magnetizations M v = (Sf) . The terms C' r and Cf, Eqs. flUEJ) are related to the 
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probabilities for the two considered spins to be either in antiparallel or in parallel Bell 
states, respectively [18]. The total amount of bipartite entanglement may be estimated 
by the so called two-tangle, namely the sum r 2 = 2 J2 r > which is related with the 
one-tangle via the monogamy inequality r 2 < T\ [19,20]. The difference T\ — r 2 is the 
so called residual tangle, while the ratio Tijrx is usually referred to as the entanglement 
ratio. 

Eqs. (J31IS1) has been originally derived [17] exploiting the symmetries of the 
Hamiltonian Eq. (JTJ), for this reason in presence of spontaneous symmetry breaking, 
i.e. for 7 > and h < 1 one has to be careful. This problem has been studied in 
Refs. [21,22]: Eqs. (EH generally hold when C" r > C' r - i.e. for h 2 + 7 2 > 1 [16] - 
while in the antiparallel region h 2 + / y 2 < 1 where C' r > C" they represent a lower bond 
for the pairwise entanglement. Moreover, from Eq. (7) of Ref. [22], one can see that 
in the asymptotic limit r — > oo they stay valid also for C' r > C", which makes the 
analysis of the long-distance concurrence, reported below, valid both inside and outside 
the factorized circle. As for the XX model, no spontaneous symmetry breaking occurs 
for whatever value of the field. 

3. Long distance concurrence 

3.1. Anisotropic case 

Let us first consider the behavior of the model for 7>0 and h > hi, in the vicinity of the 
factorized circle: we will keep fixed and finite the value of 7, and vary the field, meaning 
that we will move along horizontal lines in the h — 7 plane. 

The T = correlation functions entering the expressions of C r for the XY model 
in a transverse field are usually evaluated numerically, by computing the corresponding 
Toeplitz determinants, and cannot be written in closed form for generic r, except in the 
case of factorized ground states, where they do not depend on r. Since we are interested 
in the behavior of the concurrence as the factorized circle is approached, we fix the value 
of 7 and derive C r as a series expansion in the difference h — hi, with the factorizing 



field hi = a/1 — 7 2 . The entries of Toeplitz determinant are basically given by the well 
known G function [23] 

1 r (/i-cos0)cos(r0)+7sin0sin(r0) 
G(r,ft,7) = -/ d(p — — , 6 

^ Jo HKr,<p) 



with X(h,j; 4>) = {h — cos <fi) 2 + 7 2 sin 2 <p, and we have to expand it in the difference 
h — hi, thus obtaining 



G(r,h n ) = - 



I r* ^a/1— 7 2 — cos^j cos(r0)+7 sin <p sin(r<; 



Jo M0) 

(7) 



cos(r0) / (a/1— t 2 — cos. 
~JiW { Af(0j 
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(h-h f ) + ... 



with Af (0) = A(/if , 7; 0). The trigonometric functions of the angle rip can be represented 
in terms of powers of these functions, then, integrating by part and resumming all the 
r terms, one gets [24] 

-a r (h - h { ) + 0(h - h) 2 , 



G(r,/i, 7 ) 
G(0,h, 7 ) 



2 7 

a + ^-(h-h f ) + 0(h 

2 7 



hi) 2 



(8) 



2 ^-ft— -f 



l+2r7-(2r 2 + 1)Y 



-a 



r-2 



(h-h f )+0(h-h f 



1+7" 2 7 ( 7 +i) 2 
where the first and the third equations hold for r 7^ 0. From Eqs. ([8]), the expansion for 
the correlators are found in closed form: 



9r 



9? 



27 
1 + 7 



+ 



a 



2r+l 



2a 



(h-h f 



f -iy a 



2r-l 



fJr 



a 



2 7 



-(h 



2-, 

+ Q(h - h) 2 



+ 0(h 



h f ) 2 



„2r-l 



a a 

~ + 7 

7 7 + 1 



(h - hi) 



+ 0{h-hif 



and the magnetization along the field direction reads 

M K 



ry 1 

- + -(h-hi) + 0(h-h ( ) 2 



(9) 
(10) 

(11) 
(12) 



In the most anisotropic 7 = 1 case, it is a = and the only finite correlator up 
to the first order in (h — hi) is g® x , whose modulus gets the maximum value (i.e. 
\g* x \ = 1/4) independently of h and r; the first correction is of order (h — hi) 2 , being 
\g* x \ = l/A-(h- hf) 2 /16. Notice that Eqs. (JHU) do not hold for 7 = 0, where in fact 
they display unphysical singularities. 

Similar expressions are found for h < h{ and, by defining the distance e = \h — hf\, 
we obtain the first-order term of the expansion in e for C r , which reads 

„2r-l 



ft" 



27 



-s + 0(e 2 ) 



(13) 



and noticeably holds for whatever r. 

Let us now focus on the range R of the concurrence [16, 18], which is the distance 
between the two farthest entangled spins along the chain, i.e. 

R : C r > , Vr < R A C r = , Vr > R. (14) 

Since C r is finite for all r at the first order in e, R diverges for e — > 0. This statement 
is compatible with the behavior depicted by the exact numerical data shown in Fig. [TJ 
from the same data, we also see that for whatever r > 1, it exists a distance Eo(r, 7) 
such that C r > for < e < e (r, 7). On the other hand, since Eq. (TlBl cannot describe 
the vanishing of C r at Eo(r,'-f), in order to further analyze the behavior of R, one has 
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Figure 1. Entanglement phase diagram: The thick curve is the line where the ground 
state is exactly factorized. The regions between curves with the same drawing represent 
the areas of the h — 7 plane where C r ^ for r = 3 (dot-dashed) and r — 4 (dashed). 
C2 = only below the solid line and at h = h{ and C\ vanishes only at h = hf. The 
inset shows C r vs h for different r = 1, ...,4 (from top to bottom) at 7 = 0.5. 



to evaluate the concurrence up to the second order in e. As we are interested in the 
behavior of the long-distance concurrence, we are allowed to use the large-r asymptotic 
expressions of the correlators [23], thus finding 

^ = ^e-[A*-5A 2 (r)]e 2 + 0(e s ), (15) 

where A 2 = a 2 (7 + 3)/327 3 , and SA 2 (r) ~ 0(r~ 2 ). It is important to notice that, in 
contrast to Eq. (1131) . the above expression only holds for large r. The behavior predicted 
by Eq. (Tl5|) . though approximated, is consistent with that shown by the exact numerical 
data (see Fig. [2]). In particular, beyond the trivial zero in e = 0, Eq. (T151) has another 
zero which approximates £o( r ;7) : 

a 2r ~ l 

£o(r,7) ^4^' (16) 
where we neglected the SA 2 (r) term, which vanishes for r — > 00. Notice that the 
symmetry of Eqs. (fl5l) and ([16]) with respect to the sign of the difference h — hf arises 
from the 2 nd order expansion in h — h{ itself, that becomes more and more accurate for 
larger r (see Fig. [2]). 

For a given (large) r, Eo{r, 7) is the distance from h{ at which C r gets finite while 
approaching the factorizing field. We can rephrase this statement by saying that, for 
fixed h ^ hf, the farthest entangled spins are those whose distance r fulfills Eq. ( TT6|) . 
with 60 = \h — hf\. Therefore, if we consider r as a continuous variable, Eq. (fl6l) can 
be inverted and we obtain, for the range of the concurrence defined in Eq. (fl4l) . the 
following expression 

1 , ln(4«7A 2 ) „ , . 

R~- -lne+ \ \ ' for e^O. 17 

In cr In cr 

From the above expression we see that, for a fixed value of e, a larger anisotropy 
implies a smaller value of R. On the other hand, when the anisotropy increases, one 
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Figure 2. C r versus h — hf, for 7 = 0.5, r = 5 (top panel) and r — 8 (bottom panel): 
Comparison between the exact (full line) and approximated [Eq. (fT5j> with 5 A 2 (r) = 0] 
value (dashed line). 



should pay particular attention to the overall consistency of the reasoning, as the 
existence itself of Eo is not generally due for small r, given that Eq. ( fl7|) holds only 
for large r. Specifically, for h < hf C\ is always finite, and for h > hf both C% and 
C2 keeps finite no matter the value of the field, as seen in Fig. 1 . In particular, for 
7 = 1 the above scheme breaks down: For any finite value of the field, the only non-zero 
concurrences are those between nearest and next-nearest neighbors, as from the exact 
results by Pfeuty [25], which give 

3.2. Isotropic case 

Factorization and quantum criticality are two distinct phenomena, occurring usually 
for different values of the external magnetic field, being typically < hf < h c . When 
factorization occurs, peculiar features of the two-spin entanglement distribution are 
observed at hi where, in turn, standard magnetic observables behave quite trivially. On 
the other hand, it is just the peculiar behavior of these latter properties that signals 
the occurrence of a quantum phase transitions at h c , where two-spin entanglement 
distribution has no distinctive features. 

When factorization and quantum criticality get to coincide at hf = h c , standard 
two-point correlation functions and pairwise entanglement together signal the occurrence 
of a phenomenon which corresponds both to a factorization (though of a particular type, 
i.e. saturation) and to a quantum phase transition (though of topological type rather 
than second order). 

The specificity of the above depicted situation translates into a peculiar distribution 
of the two-spin entanglement along the chain, which gives the isotropic XX model a 
special role in the overall analysis, as shown below. 
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Let us consider Eq. (pQ) with 7 = 0, in the non-trivial quasi-ordered phase, 
h < h{ — 1: Using a procedure similar to that depicted above, we obtain the following 
closed forms for the expansions of the correlators 



9 



9r 



xx (yy) 



-1/2 



(4r 2 



L 7T 



y/2 12vr^2 



i) £ 3/ 2 + 2r(r 2 



9tt 2 



+ 0(e 



1 

4 



/2 



r 3/2 



4r 



7T 



+ — ^ 2 + 0(e 



5/2> 



(19) 



67r\/2 3tt 2 

where gr* 35 = due to the symmetry in the xy plane, and the magnetization along the 
field direction reads 



M, 



1 

— cos 

7T 



-1 



h . 



(20) 



The expansion of the concurrence consequently reads 



2^2 



.1/2 



4r 8r v / 3-(4r 2 -l)7T 



7T 



?rV3 



e 3 / 2 + 



+ 2r 



30^ + 20(r 2 -l)7r + (4r 2 -5)v / 3vr 2 ) 

457T 3 



e 2 + 0(e 



5/2N 



(21) 



From the above expressions, we see that the change of the universality class at 
7 = has drastic effect on the mechanism of rearrangement of two-spin entanglement 
along the chain. In particular, the fact that the correlation functions have an algebraic 
dependence on r, rather than the exponential one found in the anisotropic case, reflects 
in the independence of r of the first term of Eq. (|2"T|) . Thus, whatever the selection of the 
two spins in the chain they share the same amount of entanglement. Such invariance is 
quite a surprising feature and it suggests, according to the analysis proposed in Ref. [26], 
the more relevant role of multipartite entanglement in the isotropic case. 

Moreover, the comparison between Eqs. (I2T1) and (|T5l) . as well as that between the 
inset of Fig. [Hand Fig. [3j evidences an overall increase of all the concurrences {C r } in 
the most isotropic case, consistently with what is observed in the XXZ model [5]. 




C 



0.2 



0.1 



i , 1 , ^ , 1 / 11 11 1 1 

0.25 0.5 0.75 1 

h 

Figure 3. concurrence C r vs h for 7 = (XX model) and r = 1, 7 (from the 
highest to the lowest curve). 
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Let us now study R as saturation is approached from below, i.e. for h — > hj: Its 
divergence is favored by the increase of the symmetry, as testified by the singularity 
of the prefactor 1/lna 2 in Eq. ffTTl) as 7 — > 0, that signals a qualitative change in the 
behavior of R. In contrast to the anisotropic case, for 7 = we already got C r to 
order e 2 [Eq. ( |2T|) ] and we do not have to resume the 1/r asymptotic expansions of the 
correlators to evaluate the range of the concurrence. The farthest entangled spins are 
those whose distance r fulfills Eq. (I21l) =0; the latter is an equation of the 3rd order in r 
that for sufficiently small e has three real solutions, the smallest positive one is just R. 
As in previous large r studies based on the XXZ model [16, 18], we find that R diverges 
more rapidly than in the anisotropic case [Eq. (|T71) ]. namely 

R oc e- 1/2 . (22) 

For h > 1, being the ground state factorized, all the entanglement measures vanish; 
however one may fix a value of h larger than unity, and study the behavior of the 
concurrences as 7 — > 0. We therefore computed the correlators and the concurrences in 
this parameter region as a function of the anisotropy. In Fig. HJ {C r .} r . = i i ... i6 are plotted 
versus 7 at fixed field. Again we see that, while approaching the factorized ground state, 
i.e. for 7 — > 0, all the {C r } get finite and the range of the concurrence diverges with a 
logarithmic trend: R oc l/hi7. Our results show that this behavior is general for any 
h > 1 and the divergence becomes more and more pronounced for h — > 1. Eventually, 
as shown by Fig. [5j at the saturation field hf = h c = 1 the divergence of the concurrence 
range modifies its dependence on the anisotropy, being R oc I/7. Thus, the change of 
the character of the divergence of R from logarithmic to power-law observed both for 
7 = and h — >• and for h = h { and 7—^0 suggests the critical point of the isotropic 
model to represents a peculiar point in this context, as also proposed in Ref. [27], though 
in a slightly different sense. 
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4. Two-spin entanglement length £ 2SE 

The most noticeable feature of the first-order expansion in e of C r , Eq. (113j) . is the 
purely exponential dependence on r, C r ~ a 2r , which indicates that a characteristic 
length emerges in the system near factorized ground states. This length scale is 

r 1 

^ 2SE = ~ln(« 7 C r ) = ' (23) 

and we have named it two-spin entanglement length as it specifically characterizes the 
distribution of entanglement between different spin pairs along the chain. In fact, by 
looking at the expansions of the correlators, Eqs. (I MTTj) . we notice that g yy has the same 
purely exponential behavior of C r , which means that £ 2SE coincides with the standard 
correlation length along the less favored direction which, in turn, does not enter the 
characterization of the magnetic behavior. In fact, the relation between g yy and C r is 
not accidental: as a matter of fact the change of sign of g yy , namely g yy < for h < hi, 
g yy > for ft > /if, and g y r y — at /if, rules the swap between parallel and antiparallel 
entanglement, i.e. between C' r and C" [see Eqs. (HI) and (JSJ)]. So that surprisingly g yy , 
which is the less significant correlator as far as the standard magnetic properties are 
concerned, plays a relevant role in determining spin pair entanglement properties close 
to the factorization. 

The above definition Eq. fl23|) makes sense only if Eq. ffl3l) holds, i.e. for very small 
e, where the concurrence is finite for any spin pair along the chain. In other terms, when 
R diverges C r is found to decay exponentially with a characteristic length £ 2SE which gets 
larger and larger as the anisotropy decreases, finally diverging as the isotropic critical 
point (h = 1,7 = 0) is approached along a vertical line, being 

£ 2SE ~^ for 7-0. (24) 
27 

In order to further investigate the behavior of the isotropic model, let us extend 
our analysis to the vicinity of the factorized line: In the previous Section we observed 
analogies between the behavior of the two-spin entanglement close to the factorized 




Figure 5. R versus I/7 for h = h c = !• 
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^2SE 
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2 



"1 1.1 1.2 1.3 1.4 1.5 

h 

Figure 6. £ 2SE vs h at 7 = 7.5 1CP 9 . The dashed-line is the best fit f(h) = 
0.055 + 0.69/O - l) 50 . In the inset the lin-log plot of C r vs r, the slopes of the 
lines correspond to the values of £ 2SE in the main panel. 

circle (i.e. h < 1) and close to the factorized line (i.e. h > 1). In particular we found 
that also for h > 1 and 7 — > all the concurrences {C r } become finite, so that one can 
ask whether it is possible to extend the definition of £ 2SE to the region h > 1 and 7 <C 1. 
We computed C r versus r for a fixed value of 7 <C 1 and also in this case we found 
an exponential dependence on r, as shown by the inset of Fig. [61 Thus, extending the 
definition ( 1231) to the region h > 1, the two-spin entanglement length can be evaluated 
by the slope of the lines in a lin-log plot of C r vs r. As the point 7 = and h — 1 is 
approached from /i > 1, our numerical results show that 

£ 2S e oc Jj—^y for h 1+ ' ( 25 ) 

with z/ = 0.50. We notice that, by using the identity 7 = a/1 — , one may recast 
Eq. (El in the form £ 2SE oc 1/[2(1 - tf) 1 / 2 ] ~ 1/[2 V / 2(1 - h f ) 1/2 ] which reproduces the 
behavior of the above Eq. (125]) : This tells us that the way the two-spin entanglement 
length diverges while approaching the (7 = 0, h = 1) critical point does not depend on 
whether one moves /i or 7. 

The divergence of the two-spin entanglement length for 7—^0 and h — > 1 means 
that in the neighborhood of this point, not only all the concurrences {C r } are finite 
for any r, but also that the pairwise entanglement does not depend on the distance r 
between spins, as the expression for the concurrences Eq. ( |2T1) anticipated. Elsewhere 
the concurrence is either vanishing for short distances or exponentially suppressed with 
r. The behavior of £ 2SE confirms the peculiarity of the critical point of the isotropic model 
among those where the ground state of the system gets factorized. We understand such 
peculiarity as related with the fact that the critical point of the isotropic model does 
in fact coincides with saturation, i.e. with a special case of factorization. Therefore, 
while peculiar features of entanglement properties and standard magnetic properties 
are usually observed for different values of the field, hf and h c respectively, in the 
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isotropic model they occur together at the saturation field. A significant consequence 
of this feature is that, due to the lack of anisotropy in the XY plane, the XX model 
has g% x = gf y and the correlation lengths along the x and y direction are consequently 
identical: Therefore, the two spin entanglement length, that we have found to equal the 
correlation length along the y direction, in the isotropic, 7 = 0, model coincides with 
the relevant correlation length along the x direction, with which consequently shares 
the divergence at h c = h { . 



5. Residual entanglement and entanglement ratio 

The amount of entanglement stored between two spins far apart in the chain does not 
only rely on the distribution of concurrences {C r }, but also on the total entanglement of 
the system and on whether it is bipartite or multipartite. A simple way to investigate this 
issue is to evaluate the one-tangle T\, the residual tangle T\ — r 2 , and the relative weight 
of the pairwise entanglement through the entanglement ratio 72/71, in the neighborhood 
of the factorized circle. For < 7 < 1, using the definitions of one- and two-tangle and 
the expressions for the magnetizations and concurrences, one obtains 

(l_ 7 )( 3 + 7 ) 



*7 3 (l + 7) 



e* + 0(e 6 ), (26) 



(l_ 7 )2( 2 + 7 ) 

ri ~ T2= 87^(1+7) " + ° (£h (27) 

^ = ^- + 0(e). (28) 
ri 3 + 7 

For 7 = 0, as noticed above, the correlators do not decay exponentially with r, 
but they rather follow a power law, as a consequence of the quasi-long range order 
characterizing the ground state for h < 1 [compare Eqs. (!9lfTT!) with Eqs. (fl8|) and (|T9l) ]. 
This behavior, which reflects on the concurrences Eq. ( J2T1) . makes it cumbersome to 
evaluate the sum in the two-tangle expression. In fact, it is not difficult to show that 
t 2 oc e 1 / 2 , but in order to obtain the proportionality coefficient with good accuracy one 
should retain several terms in the small e expansion of C r . For this reason we preferred 
to numerically compute r 2 and t 2 /tx close to the factorizing field up to e = 10~ 6 and 
we verified that the first term of Eq. ( 1281) holds also for 7 = 0. 

From the analysis of the above expressions, we notice that the larger the anisotropy 
the higher the relative weight of the pairwise entanglement close to the factorized ground 
state. In particular, in the pure Ising limit 7 = 1 the entanglement in the ground state is 
totally stored in pairwise form up to order e 2 , but both the total entanglement and the 
residual tangle are strongly suppressed, being n = e 4 /32 + 0(e 6 ) and T\ — r 2 = e 6 /64. 
In the opposite limit, small values of the anisotropy 7 favor the presence of multipartite 
entanglement and reduce the relative weight of the two-spin entanglement. In order 
to reconcile this last statement with the fact that, as shown in Section 13. 2\ a smaller 
anisotropy implies larger C r , one should notice that, as the anisotropy of the model 
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decreases, the one-tangle becomes larger and larger and, for 7 = it is 

n = 1 - 4M 2 2 

= ^VJ.L + 0(^) , (29) 

7T 7T Z 

where M z is given by Eq. (120]) . Thus, even if in the 7 — > limit the fraction of pairwise 
entanglement reduces, the two-spin entanglement takes advantage of the overall increase 
of the total amount of entanglement stored in the ground state. 

6. Conclusions 

In this paper we studied how the two-spin entanglement is distributed along the 
spin chain described by the Hamiltonian (CD), focusing our attention on the possible 
connections between the pairwise entanglement spreading and the symmetry of the 
model. In particular we have considered the neighborhood of factorized ground states, 
where all the concurrences get finite, no matter the distance between the two spins, 
and the range of the concurrence diverges [see Eqs. ([IT]) and (122]) ]. We have derived 
closed analytical expressions for correlators and concurrences in the neighborhood of 
the factorized circle: Using these formulas we have derived an analytical expressions for 
the range of the concurrence R, and shown that, whenever R diverges, a characteristic 
length scale £ 2SE , Eq. (123]) naturally emerges in the system. This two-spin entanglement 
length, which defines the pure exponential decay of C r , evidenced by Eq. tjTB"]) . is finite 
for finite anisotropy and diverges as the isotropic critical point is approached, where it 
in fact coincides with the magnetic correlation length on the plane perpendicular to the 
applied field. 

Our description shows how the two-spin entanglement distribution along the chain 
evolves, while moving from the Ising (7 = 1) to the isotropic (7 = 0) model: In fact, for 
7 — ► 1, it results that a finite amount of entanglement can be stored between nearest 
neighbor spins only; in addition in the pure Ising case the range of the concurrence is 
always finite. On the other hand, the fully isotropic exchange interaction of the 7 = 
model evidently favors the entanglement rearrangement also between distant spins. 

Remarkably enough, despite the concurrence is either vanishing for short distances 
or exponentially suppressed with r in any point of the parameter plane h — 7, in the 
isotropic case, close to h — h c — 1 all the concurrences {C r } are finite and their value is 
independent of the specific pair of spins along the chain whatever the distance r between 
them. This fact together with the change of the divergence character of the concurrence 
range testifies the special role played by the critical point of the isotropic model in the 
distribution of the entanglement between two spins of the system. 
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